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ABSTRACT 

Context. A small fraction of upper main sequence stars have strong, highly structured magnetic fields. The origin and evolution of 
these fields are not adequately understood. 

Aims. We are carrying out a survey of magnetic fields in Ap stars in open clusters in order to obtain the first sample of magnetic upper 

main sequence stars with precisely known ages. These data will constrain theories of field evolution in these stars. 

Methods. A survey of candidate open cluster magnetic Ap stars was carried out using the new ESPaDOnS spectropolarimeter at the 

CFHT This instrument provides an alternative to the FORS 1 spectropolarimeter used up to now for this survey. 

Results. We have obtained 44 measurements of the mean longitudinal fields <S ; > of 23 B6 - A2 stars that have been identified as 

possible Ap stars and that are possible members of open clusters, with a median uncertainty of about 45 G. Of these stars, 10 have 

definite field detections. Nine stars of our sample are found not to be magnetic Ap stars. These observations significantly increase the 

information available about low-mass stars near the TAMS compared to our previous sample. 

Conclusions. We find that ESPaDOnS provides field measurements comparable to those that we have previously obtained with 
FORS 1 , and that these data also contain a large amount of useful information not readily obtained from lower resolution spectropo- 
larimetry. With the new data we are able to expand the available data on low-mass, relatively evolved Ap stars, and identify more 
robustly which observed stars are actually magnetic Ap stars and cluster members. Re-analysis of the enlarged data set of cluster Ap 
stars indicates that such stars with masses in the range of 2 - 5 M Q show RMS fields larger than about 1 kG only when they are near 
the ZAMS. The time scale on which these large fields disappear varies strongly with mass, ranging from about 250 Myr for stars of 2 
- 3 M Q to 15 Myr for stars of 4 - 5 M Q . Our data are consistent either with emergent flux conservation for most (but not all) Ap stars, 
or with modest decline in flux with age. 

Key words. Stars: magnetic fields - Stars: chemically peculiar - Stars: evolution - Polarization - Techniques: polarimetric 



1. Introduction 

The discovery that some A- and B-type stars of the upper main 
sequence possess strong (kG) magnetic fields was made more 
than 60 years ago (Babcock I1947I ). Since that time, there has 
been much investigation into the occurrence and structure of 
these fields, but much about them remains unknown. In partic- 
ular, the evolution of these fields with time is a major area of 
ignorance, and until quite recently, almost no observational con- 
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straints have been available to guide and test theoretical ideas 
about field evolution. 

In the past few years the observational situation has changed 
dramatically with the development of powerful new spectropo- 
larimeters such as the one developed for the MuSiCoS spectro- 
graph at the TBL, and FORS1 on the ESO VLT. With FORS1 
one can obtain routine magnetic measurements of main sequence 
stars fainter than my ~ 10 m , and these new instruments have also 
made possible a general decrease in field measurement uncer- 
tainty by a factor of several. 

The new generation of spectropolarimeters has for the first 
time made it practical to study observationally the evolution of 
magnetic fields in magnetic Ap stars with time, by exploring 
a large sample of such stars that are members of open clus- 
ters. Cluster studies open the possibility of obtaining a sample 
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of magnetic stars with reasonably accurate ages, thus providing 
a statistical time sequence of upper main sequence magnetism. 
Photometric and spectroscopic surveys have identified a substan- 
tial sample of probable Ap stars in such clusters which can be 
surveyed for fields with the new instruments, and the member- 
ship of these stars in their clusters can be established much more 
securely than in the past thanks to recent advances in astrometry, 
especially the Hipparcos and Tycho-2 projects. 

These new possibilities have stimulated a major survey of 
magnetic fields in cluster Ap stars, whose first results have been 
reported by Bagnulo, Landstreet, Mason et al. (2006 hereafter 
Paper I) and analysed by Landstreet, Bagnulo, Andretta et al. 
(2007 hereafter Paper II). In the work completed so far, we 
have obtained and collected magnetic measurements of 8 1 stars 
which are both probable magnetic Ap stars and probable clus- 
ter members. The absolute (log) ages of the stars in this sample 
are known with accuracy in the range of 0.05 - 0.2 dex (the age 
uncertainty of the clusters themselves). The fraction of the main 
sequence lifetime elapsed (the "fractional age") of each of the 
stars of this sample is known with a precision that is usually less 
than half the fractional age itself, and thus for stars of small frac- 
tional age, the fractional age is also quite precise. 

From this sample, we have established conclusively that 
fields are present in Ap stars more massive than 2M from the 
ZAMS onwards. We have also shown that stars more massive 
than about 3M have a substantial decline in field strength after 
the age of about 3 10 7 yr, while less massive stars show no strong 
evidence of change of field strength even as long as 10 9 yr after 
reaching the ZAMS. 

Recently the powerful new spectropolarimeter ESPaDOnS 
has become available at the CFHT This instrument combines a 
large resolving power with high throughput and wide spectral 
coverage, and is the most powerful instrument for field measure- 
ment in stars with reasonably rich spectra and low v e sin i that 
is currently available. We have used ESPaDOnS to extend the 
survey described and analysed in Papers I and II. 

This extension of our previous survey has several aims: (1) to 
expand the survey into the northern celestial hemisphere, adding 
clusters not visible from ESO; (2) to increase the number of ob- 
servations of stars for which only one or two magnetic measure- 
ments are available, thereby improving the precision of B ims val- 
ues available for stars of the survey; (3) to search for fields in 
probable magnetic Ap stars for which no fields have yet been 
detected by the available observations, exploiting the superior 
ability of ESPaDOnS to detect weak fields, and to detect fields 
even when the mean longitudinal field (B z ) is close to zero (see 
discussion below); and (4) to compare the field measurement ac- 
curacy obtained with high- and low-resolution spectropolarime- 
try for various kinds of magnetic Ap stars. 

The next section of this paper discusses the observational 
strategy employed, the data reductions, and the choice of stars to 
observe. Sect. 3 describes the transformation of the polarisation 
measurements into measurements of (B z ), and tests of the data 
quality. Sect. 4 presents the new field measurements obtained. 
In this section we also compare the observational efficiencies of 
various instruments. Sect. 5 contains a star-by-star assessment of 
the new information obtained from the ESPaDOnS observations. 
In Sect. 6 we discuss general features of the new addition to our 
data set of cluster magnetic Ap stars. In Sect. 7 the new cluster 
field measurements are integrated into our larger survey, and the 
conclusions of Paper II are reviewed in the light of our new data. 
Finally Sect. 8 summarises the conclusions of the paper. 



2. Instrumentation, observations and data 
reduction 

2.1. ESPaDOnS 

The instrument used for this portion of our survey of mag- 
netic Ap stars in clusters is ESPaDOnS, the new cross- 
dispersed echelle spectropolarimeter built for the Canada- 
France-Hawaii telescope. The instrument is conceptually similar 
to the MuSiCoS spectropolarimeter which has been extensively 
used for high-precision magnetic measurements (e.g. Wade et 
al. [2000a 2000b), but has a factor of about 30 times higher effi- 
ciency. 

The polarisation analyser is located at the Cassegrain focus 
of the telescope. The stellar image is formed on an aperture fol- 
lowed by a collimating lens. The beam then passes through a ro- 
tatable A/2 waveplate, a fixed A/4 waveplate, a second rotatable 
A/2 waveplate, and finally a small-angle Wollaston prism, fol- 
lowed by a lens which refocuses the (now double) star image on 
the input of two optical fibres. This relatively complex polarisa- 
tion analyser is necessary because one of the fundamental design 
parameters for ESPaDOnS was very wide wavelength coverage 
(approximately 3700 A to 1.04 jj.m). To have waveplates which 
are approximately achromatic over this wide range, ESPaDOnS 
uses Fresnel rhombs. A single Fresnel rhomb acts as a A/4 plate, 
but deviates the beam, while two Fresnel rhombs in series form 
a A/2 plate without beam deviation. To minimise mechanical 
complications, only the double (non-deviating) Fresnel rhombs 
are allowed to rotate; the configuration chosen is the minimum 
which allows one to analyse all of the Stokes polarisation com- 
ponents (Q, U, V) by appropriate orientation of the axes of the 
successive waveplates. 

The two output beams from the Wollaston prism, which have 
been analysed into the two components of linear or circular po- 
larisation as desired by appropriate waveplate settings, are then 
carried by the pair of optical fibres to a stationary and thermally 
buffered (but not yet temperature-controlled) cross-dispersed 
spectrograph where two interleaved spectra are formed, covering 
virtually the entire desired wavelength range with a resolving 
power of R « 65 000. The / component of the stellar Stokes 
vector is formed by adding the two corresponding spectra, 
while the desired polarisation component (Q, U, V) is obtained 
essentially from the difference of the two spectra. To minimise 
the systematic errors due to small misalignments, differences in 
transmission, effects of seeing, etc., one complete polarisation 
observation of a star consists of four successive spectra taken 
with four different positions of the two rotatable Fresnel 
rhombs; for the second and third exposures, the waveplate 
settings exchange the right and left circularly analysed beams 
with respect to the handedness measured in the first and fourth 
exposures (cf. Donati et al. 1997; Donati et al. in preparation; 
http:www.cfht.hawaii.edu/Instruments/Spectroscopy/Espadons/) 

2.2. Observations 

The stars observed were selected from a large data-base we have 
assembled (see Papers I and II) of possible magnetic Ap stars 
that are also possible members of open clusters within about 1 
kpc of the Sun. As discussed in Paper II, we have used avail- 
able astrometric and photometric data, especially those from the 
Hipparcos (ES A fTWTj l and Tycho-2 (H0g et al. |2000al [2000bb 
catalogues, to clarify as far as possible which magnetic Ap stars 
are probable or definite cluster members. Here we use essentially 
the same criteria for cluster membership that were discussed in 



J.D. Landstreet, J. Silaj, V. Andretta et al.: A survey of open cluster Ap stars with ESPaDOnS 



3 



Paper II. We selected stars from the data-base that were observ- 
able in July 2005, when these observations were carried out. 
In part, we observed stars that are not observable from ESO, 
where we have obtained the previous observations for this sur- 
vey. However, a number of observations were obtained of stars 
for which previous measurements are available, either by us or 
by others. 

As discussed in detail in Paper II, the classification of most 
of the stars observed as magnetic Ap stars is still fairly uncertain, 
because identification of cluster Ap stars is usually based only on 
narrow-band photometry and/or low dispersion, low signal-to- 
noise spectra. Thus we expected that some of the stars observed 
would turn out not to be magnetic Ap stars at all. 

Our previous observations for this survey were carried 
out using the low -resolution spectropolarimeter FORS1 at the 
ESO VLT. With R ~ 10 3 , FORS1 magnetic field measure- 
ments give high weight to the circular polarisation signal in the 
Balmer lines, which are resolved even at this low dispersion. In 
cases where a relatively rich metallic spectrum is present, the 
FORS1 measurements can be made significantly more precise 
by exploiting the polarisation signal present in the (barely re- 
solved) metallic line spectrum, but to first order the precision 
of FORS 1 magnetic field measurements is not sensitive to the 
details of the metallic line spectrum. In contrast, measurements 
with ESPaDOnS, like those obtained with MuSiCoS, have high 
enough resolving power to fully exploit the information content 
of the metallic line spectrum. ESPaDOnS field measurements 
have standard errors which decrease strongly with decreasing 
v e sin i and with increasing richness and strength of the metallic 
line spectrum (cf. Landstreet 1982; Shorlin et al. 2002). To fully 
exploit this dependence, and thus to obtain the most precise mea- 
surements possible, we have preferentially observed relatively 
cool magnetic Ap stars, and have selected our targets for low 
v e sin i where available data allowed us to select targets with this 
criterion. 

The list of stars observed is given in Table Q] which contains 
the name of each star, the cluster to which it may belong, our 
assessment of the probability that the star is actually a cluster 
member (Y = definite member, P = probable member, ? = ques- 
tionable, assessed as in Paper II), the spectral type based on a 
combination of published types with the results of our own ex- 
amination of individual / spectra, and the v e sin i value deduced 
from our ESPaDOnS spectra (the uncertainty in these values is 
approximately the larger of ±10% or ±2 km s ). 

2.3. Data Reduction 

The observed polarised spectra were reduced to ID using the 
dedicated software LibreEsprit, provided by J.-F. Donati for 
treatment of ESPaDOnS data. LibreEsprit subtracts bias, locates 
the various spectral orders on the CCD image, measures the 
shape of each order and models the (varying) slit geometry, 
identifies comparison lines for each order and computes a global 
wavelength model of all orders, performs an optimal extraction 
of each order, and combines the resulting spectra (in groups 
of four, corresponding to the four sub-observations with the 
four different Fresnel rhomb settings) to obtain intensity (/) and 
circular polarisation (V) spectra. The V spectrum normally has 
the continuum polarisation removed, as this arises mainly from 
instrumental effects and carries little information about the star. 
Each spectrum is corrected to the heliocentric frame of refer- 
ence, and may optionally be divided by a fiat field and be approx- 
imately normalised (see Donati et al. 1997 and Web pages at 
http://www.cfht.hawaii.edu/Instruments/Specti-oscopy/Espadons) 



Table 1. Stars observed with ESPaDOnS for magnetic fields 



Star 


Cluster 


Memb 


Spectrum 


v sin i 
(kms" 1 ) 


HD 16605 


NGC 1039 


P 


Alp SiCrSr 


18 


HD 16728 


NGC 1039 


Y 


B9V 


75 


HD 19805 


a Per 


Y 


B9.5V 


8 


HD 108945 


Coma Ber 


Y 


A2p SrCr 


65 


HD 144661 


Upper Sco 


Y 


B8p He-wk 


45 


HD 153948 


NGC 6281 


P 


Ap Si 


80 


HD 317857 


NGC 6383 


? 


AlIVp 


26 


HD 318100 


NGC 6405 


Y 


B9p 


42 


HD 320764 


NGC 6475 


P 


AlVan 


225 


HD 162576 


NGC 6475 


Y 


B9.5p SiCr 


30 


HD 162588 


NGC 6475 


Y 


ApCr 


73 


HD 162630 


NGC 6475 


P 


B9V 


51 


HD 162656 


NGC 6475 


Y 


B9V 


45 


HD 162725 


NGC 6475 


Y 


AOp SiCr 


32 


HD 169842 


NGC 6633 


Y 


Alp SrCr 


45 


HD 169959A 


NGC 6633 


? 


AOp Si 


57 


HD 170054 


NGC 6633 


Y 


B6IV 


26 


HD 170860 


IC 4725 


P 


B9IV/Vp 


85 


HD 172271 


IC 4756 


P 


AOpCr 


107 


HD 205073 


NGC 7092 


Y 


Al 


16 


HD 205331 


NGC 7092 


Y 


A1V 


51 


BD+49 3789 


NGC 7243 


P 


B7p Si 


81 


HIP 109911 


NGC 7243 


Y 


AOVp Si 


59 



A check spectrum called the N spectrum, computed by 
combining the four observations of polarisation in such a way 
as to have real polarisation cancel out, is also calculated by 
LibreEsprit. The N spectrum tests the system for spurious po- 
larisation signals. In all of our observations, the N spectrum is 
quite featureless, as expected. 

The final spectra consist of ASCII files tabulating /, V, N, 
and estimated uncertainty per pixel as a function of wavelength, 
order by order. We have kept both normalised and un-normalised 
versions of each spectrum. 

3. Measurement of magnetic field strength 

Two methods allow us to detect magnetic fields in the stars ob- 
served. First, we can use the V spectra to measure the mean lon- 
gitudinal magnetic field strength (B z ) of each star at the time of 
observation. This is the conventional measure of field strength 
normally used for detection of fields in magnetic Ap stars (e.g. 
Landstreet 1982 ). However, because of the high value of resolv- 
ing power, we can also examine spectral lines for the presence 
of circular polarisation signatures; Zeeman splitting combined 
with Doppler broadening of lines by rotation can lead to non- 
zero values of V within spectral lines even when the value of 
{B z ) is close to zero. This possibility substantially increases the 
sensitivity of our measurements as a discriminant of whether a 
star is in fact a magnetic star or not, as discussed by Wade et al. 
d2000al l and Shorlin et al. (120021 . 

Because ESPaDOnS is a high-resolution instrument and be- 
cause the stars observed are relatively faint (typical magnitudes 
my between 6 and 10), the signal-to-noise ratio (SNR) in indi- 
vidual spectral lines is not very high, and field measurements 
made with single lines are not very precise. In this situation, 
we employed the technique of Least Squares Deconvolution 
(LSD), developed into a practical strategy by Donati et al. d 19971 ) 
and employed extensively in the analysis of magnetic observa- 
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tions made with the MuSiCoS spectropolarimeter (cf Wade et al. 
I2000ab . 

In the LSD technique, we start with a list of spectral lines 
expected in a normal or a magnetic Ap star of a temperature 
near that of the star observed. This line list is obtained from 
the Vienna Atomic Line Database (VALD; Piskunov et al. 1995 ; 
Ryabchikova et al. 1999; Kupka et al. 1999), and is selected by 
specifying a suitable chemical composition, wavelength range, 
and minimum line depth. For the stars in our sample, two dif- 
ferent abundance tables were employed. An Ap abundance ta- 
ble was created by assigning the iron-peak elements abundances 
lOx solar abundance, except for chromium which was increased 
to lOOx solar. A solar abundance table is the default provided 
by VALD when a request is made, and was used to make the re- 
maining lists. VALD line lists were requested to include all lines 
in the ESPaDOnS spectral range (3700-10000 A) with depths 
greater than 10% of the continuum (in a synthetic spectrum with 
v e sin / = 0). 

Fundamental parameters (effective temperature T e $ and sur- 
face gravity) for the stars observed were determined primarily 
from Stromgren and/or Geneva photometry, as described in de- 
tail in Paper II. Neither Stromgren photometry nor Geneva pho- 
tometry is available for the stars HD 169842, HD 172271 and 
HD 317857, so the effective temperatures were estimated from 
Johnson UBV photometry or the observed spectra. Initial selec- 
tion of an appropriate LSD line list for each star was made based 
on these values. Experimentation showed that the highest SNR in 
the measured field (B z ) generally resulted from the use of a line 
list selected with a temperature close to, or up to about 10 - 15% 
hotter than, the estimated T e ff of the star. (Using a line list ap- 
propriate for a slightly hotter star seems to eliminate a number 
of weak lines that contribute more to the noise than to the signal.) 

For each spectral line in the line list appropriate to a par- 
ticular star, a small window is cut out of the observed / and V 
spectra centred on the laboratory wavelength of the line, and the 
local wavelength scale is converted to a velocity scale relative to 
the nominal central wavelength. Conceptually, one may imag- 
ine that all the small / and V spectral windows are then aver- 
aged. (In practice, model intensity and polarisation signatures 
are fit to the ensemble of single line data windows, rather than 
performing direct averaging.) In this average, each window is 
given a weight proportional to the line depth, the wavelength, 
and the Lande factor of the line. Since the / and V signatures 
in individual spectral lines are very similar to one another, and 
scale essentially with the weighting factors used in averaging, 
the resulting mean intensity and polarisation line profiles may 
be analysed as if they represented a single real spectral line hav- 
ing a much higher SNR than the individual observed lines. The 
improvement in SNR from single lines to the LSD line is by 
roughly the square root of the number of lines used, so the SNR 
in the LSD mean line is typically of order 30 times higher than 
in single lines. 

The LSD line is then analysed in two ways. First, the value 
of (B-) is computed using the normal expression (e.g. Donati et 
al. [19971 

f vV(v)dv 

<A> = -2.1410 12 J - , (1) 

Azc J [I c - I(v)]dv 

where (B z ) is in G, z is the mean Lande factor of the LSD line 
(typically about 1 .25), and A is the weighted mean wavelength 
of the LSD line in A (typically about 5200 A). The limits of in- 
tegration were chosen visually for each star to coincide with the 
apparent edges of the LSD / and V lines; using a smaller window 




0.6- CrII 



4560 4570 4580 4590 

Wavelength (A) 

Fig. 1. A sample region of the / spectrum of HD 169842 with 
identification of strong lines (lower solid line) is compared to 
a synthetic spectrum (dashed line), for which Cr has been en- 
hanced by a factor of about 80. Above is the same region of the 
V spectrum, multiplied by 5 and displaced upward by 1.15 for 
clarity. 



would neglect some of the signal coming from the limb of the 
star, while a window larger than the actual line would increase 
the noise without adding any further signal, thus degrading the 
SNR below the optimum value achievable. If the resulting value 
of (B z ) is significantly non-zero (say by 4 - 5cr B in a single ob- 
servation, or by 3crg in two or more observations), we conclude 
that a magnetic field has been detected. 

In addition, the LSD line profile itself is examined. If the 
value of V rises to a statistically significant value inside the line, 
while always remaining insignificant in the neighboring contin- 
uum, we conclude that a field is detected. As mentioned above, 
this may occur even if (B z ) is not significantly different from 
zero. On the other hand, from our experience with searches for 
weak fields in known Ap stars (cf Auriere et al. 120071 1. a star 
with a favourable spectrum (hundreds of lines and v e sin i less 
than about 50 or 60 kms~') which shows no significantly non- 
zero values of V within the LSD spectral line in two indepen- 
dent high SNR observations is unlikely to host a magnetic field. 
Furthermore, since observations with ESPaDOnS yield high- 
resolution / spectra (as well as V spectra), we may test the sus- 
picion that a star is in fact not a magnetic Ap star at all by ex- 
amining the / spectrum. The combination of no significant V 
signature with a normal / spectrum has confirmed that several of 
the stars selected for observation are indeed not Ap stars. 

As an example, small segments of the ESPaDOnS / and 
V spectra of HD 169842 from JD 2453566.902 are shown in 
Figure Q] The / spectrum is compared to a synthetic spectrum 
computed using the magnetic line synthesis programme Zeeman 
(cf. Landstreet 1988) with (nearly) solar abundances except for 
Cr, which is enhanced by 1 .9 dex. Note that no trace of Zeeman 
polarisation is seen in individual spectral lines; in this star the 
Zeeman signal is too small relative to the noise level in V to be 
directly visible. The LSD / and V spectra computed from these 
data are shown in Figure [2] Averaging over 4836 spectral lines 
has reduced the noise in the LSD V spectrum by a factor of or- 
der 30 relative to the noise level in individual spectral lines, al- 
lowing the magnetic field signature in V to be clearly seen with 
an amplitude well above the small noise level outside the line. 
(The number of lines used in the LSD mask is a strong func- 
tion of both r e ff and degree of peculiarity. The Ap masks used 
have about 8300 lines for a star of ~ 8500 K, decreasing 
smoothly to about 1600 lines at ~ 15000 K. Typically the 
masks for normal composition have about 1/3 as many lines as 
the Ap masks.) 
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Fig. 2. The LSD / (lower) and V (upper; multiplied by 25 and 
displaced upwards by 1.02 for visibility) spectra of HD 169842 
from the same night as Figure Q] 
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3.1. Verification 

One important kind of verification of the performance of 
ESPaDOnS is to confirm that it does not produce spurious mag- 
netic field detections. We have tested this in two ways. First, we 
have made a precise field measurement (-9 ± 19 G) of the ap- 
parently non-magnetic star a And (= HD 358), which, in agree- 
ment with numerous other measurements made with a variety 
of instruments, shows no significant field (Wade et al. 2006). 
Secondly, we find in the present data a number of stars which 
turn out (on the basis of our new high-resolution / spectra) to be 
normal rather than magnetic Ap stars; invariably such stars also 
do not show any significant magnetic field. 

We have also verified that our field measurements are consis- 
tent with earlier work by observing the longitudinal field of the 
well-studied magnetic Ap star a 2 CVn during three nights. These 
data are presented in Table [2] We have computed phases (at the 
mid-point of each observation) according to the ephemeris of 
Wade et al. (2000b), and we compare the three new measure- 
ments with their magnetic curve in Figure [3] It is clear from the 
figure that the new observations are fully consistent with the pre- 
vious (MuSiCoS) data obtained and reduced in a very similar 
way, but the ESPaDOnS data are about a factor of five more pre- 
cise. 



Table 2. Observations of standard star a 2 CVn 



HJD 


SNR 




phase 


(B z )i 


--o- B 












(G) 


2453569.748 


680 


0.570 


±0.01 


589 : 


t 11 


2453570.747 


1040 


0.752 


±0.01 


-619: 


t 10 


2453571.738 


580 


0.934 


±0.01 


-739± 11 



4. Results 

The field measurements obtained during the ESPaDOnS survey 
are presented in Table [3] This table lists the stars for which field 
measurement was successful; the heliocentric Julian date (HJD) 
of each observation; the signal-to-noise ratio (SNR) achieved per 
A of spectrum, derived from the value reported in the output file 
of LibreEsprit for the N spectrum of order 45 at 5030 A; the 
measured value of the stellar heliocentric radial velocity v r (in 
km s~' ; generally accurate to about ±1-2 km s ); the measured 
mean longitudinal field strength (B z ) (in G) with its standard er- 
ror (computed by treating each integration pixel as independent, 
using the uncertainty per pixel of V assigned by the LSD pro- 
cess); and finally a flag (fid) that indicates whether a field has 



-1000 1 1 1 1 1 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

Phase 

Fig. 3. Phased magnetic mean longitudinal field variations of a 2 
CVn as reported by Wade et al. [2000b (open circles). Our three 
new data points (filled circles) are those with very small error 
bars. 



been definitely detected (DD) or not (blank). Note that the ra- 
dial velocity has been measured from the LSD spectrum by bi- 
secting the line wings close to where they reach the continuum. 
The precision of such a measurement is limited primarily by the 
width and shallowness of the line profile and occasionally by the 
unusual shape; the absolute precision of ESPaDOnS radial ve- 
locities for very sharp-line stars is a few hundred ms _1 , and the 
stability is on the order of a few tens ofms" 1 . 

As discussed above, a field may be detected through the V 
signature even if (B z ) is close to zero, as is the case for ex- 
ample of one observation each of HD 108945, HD 162576 and 
HD 162725. On the other hand, the field of a star which is actu- 
ally a magnetic Ap star may not be detected, for example if the 
field is particularly weak (of order 100 G or less), if the spectrum 
is strongly broadened by rotation, if r er j is high and there are rel- 
atively few metal lines, or if the SNR of the spectrum obtained 
is too low. 

It is of interest to compare the relative measurement efficien- 
cies of magnetic measurements made with low and high spectral 
resolution. We can obtain such a comparison by examining mea- 
surements of the same star made with FORS1 and ESPaDOnS. 
Fortunately we have in our data sets several stars observed with 
both instruments. 

In Table [4] we compare measurements of the same star made 
with FORS1 and with ESPaDOnS. This table lists individual 
stars, magnitude my, the important parameters r er y (determined 
as in Paper II) and v e sin i which have a strong influence on the 
precision of ESPaDOnS measurements, and then for specific 
measurements with FORS1 or ESPaDOnS we list HJD of ob- 
servation, the total time required for an observation (including 
telescope setting time) At, the resulting field measurement with 
its uncertainty, and whether a field was detected (for FORS1, 
as in Paper I, nnn means no detection, DnD means definite de- 
tection only in Balmer lines, and DDD means detection in both 
Balmer lines and the metal spectrum; for ESPaDOnS measure- 
ments, ND means no detection, DD means definite detection). 

Note that overheads for a single magnetic field measure- 
ment with FORS1 total about 20 min (in "fast" mode). For stars 
brighter than my = 8 or 9, the integration time is small com- 
pared to the overheads, so the time-on-target for such bright stars 
will be independent of magnitude. (The justification for using 
FORS 1 in this inefficient mode was its uniqueness in the south- 
ern hemisphere. FORS1 showed its full capabilities in the de- 
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tection of a field in a broad-lined Ap star of my - 12.9; see 
Bagnulo et al. 2004 ) The overheads for a full ESPaDOnS mag- 
netic measurement total about 9 min, so the time-on-target for 
ESPaDOnS exposures is dominated by the actual exposure time 
for stars fainter than about nty = 6. 

In examining Table @] it is worth recalling that because 
low-resolution spectropolarimetry, such as the FORS 1 measure- 
ments, relies heavily on the polarisation in the wings of the 
Balmer lines, which do not vary greatly over the spectral range 
of our survey, the measurement uncertainties obtained using 
this class of techniques tend to fall within a moderate range. 
Inspection of Table A. 3 of Paper I shows that most of the mea- 
surements there have standard errors between about 40 and 90 G, 
while the most accurate measurements have cr B ~ 30 G.This is 
also true of the subset of data in Table|4] where all the uncertain- 
ties are in the range of 33 - 77 G. In contrast, the uncertainties 
of high-resolution spectropolarimetry, such as the ESPaDOnS 
data, are significantly more scattered. The precision of a high- 
resolution Zeeman measurement depends strongly on the num- 
ber and strength of spectral lines, on v e sin z, and on the total 
signal collected. Very precise measurements (erg as low as a few 
G) are possible for stars having many sharp and deep lines (com- 
pare Wade et al. l2006t . More than 1/3 of the ESPaDOnS mea- 
surements have ctb < 30 G, while about 20% of our observations 
have <tb > 100 G; the largest <x B s in Table[3]are over 400 G. 

Examination of Table [4] reveals that in similar Afs, 
ESPaDOnS and FORS1 often obtain measurements with 
roughly similar field strength uncertainties. The ESPaDOnS 
measurements are substantially less precise primarily for stars 
with values of v e sin; above 60 or 70 kms~ (HD 153948), for 
stars with T eff above perhaps 12000 K (HD 170054), or for 
stars for which insufficient signal was obtained due to faint- 
ness, weather or other conditions. On the other hand, in cases 
of particularly favourable circumstances (cool stars with low 
v e sin z), the <x B s can be very low, even below 20 G (HD 19805, 
HD 162725, HD 205073). 

The reason that an instrument such as ESPaDOnS, on a 3.6- 
m telescope, can be competitive with one on an 8-m telescope 
is partly a consequence of the way in which the high-resolution 
instrument can exploit more fully the information content of the 
metallic-line spectrum of the star observed. However, other fac- 
tors are also involved. In the particular case of FORS1, obser- 
vations of bright stars (my < 10) have a very low duty cycle. 
Most of the observing time is spent setting on target and clos- 
ing the loop on the adaptive optics, and on reading out the CCD. 
The actual shutter time on such a star is only a few minutes out 
of the roughly 30 min needed for a typical bright star obser- 
vation. For stars fainter than 10 m , for which longer integration 
times are required (e.g. 2 hr for the my ~ 13™ star NGC 2244- 
334; Bagnulo et al. 2004), the instrument on the larger telescope 
achieves a much higher duty cycle, and the aperture advantage 
is much more important. 



5. Individual stars 

Because each ESPaDOnS observation provides not only a mag- 
netic measurement but also a high-dispersion spectrum covering 
some 7000 A through the visible and near-IR, these data can 
provide a number of useful pieces of information. We obtain a 
measurement of the mean longitudinal field as well as supple- 
mentary information about the presence or absence of a mag- 
netic field from the behaviour of V within the LSD / line profile. 
We can also measure the projected rotation velocity v e sin i, and 



the stellar radial velocity v r . The value of v e sin i provides a con- 
straint on possible rotation periods, and v r provides a potentially 
useful test of cluster membership. We can use the shape of the 
LSD / profile, and its possible variations if we have multiple ob- 
servations, as valuable additional clues to test whether a doubtful 
star is actually a magnetic Ap star. Finally, we can examine (and 
even model) the observed / spectrum to characterise the chemi- 
cal composition of the atmosphere in order to securely establish 
peculiarity or normality. In this section we discuss the informa- 
tion obtained about each of the stars studied. 

5.1. HD 16605 = NGC 1039-154, A1p SiSrCr 

The large field of this star was first reported by Kudryavtsev 
et al. d2006l l. Our two measurements confirm their result. All 
five available measurements show negative fields, ranging from 
-840 to -2430 G. The observed v e sin /= 18 km s -1 (see Table[Q) 
suggests a rotation period of a few days. The star is certainly a 
magnetic Ap star. 

Proper motion data (see Paper II) for this star indicate proba- 
ble membership in NGC 1039, as does the star's position near 
the cluster isochrone in the HR Diagram. The radial velocity 
measured for this star is about -1 kms~'. Our two observations 
show no variation of this quantity, and the three spectra obtained 
by Kudryavtsev et al. (2006) yield v r values that are consistent 
with our measurements (Kudryavtsev, private communication). 
Thus it appears that this star is probably not a spectroscopic bi- 
nary. However, the radial velocity of the cluster is reported by 
Kharchenko et al. (2005) to be -15 + 2.3 kmr 1 . Furthermore, 
nearly 50 radial velocities have been measured in low-mass main 
sequence stars in this cluster by Jones et al. (1997); omitting 
two stars with strongly discordant radial velocities, the median 
cluster velocity is about - 1 km s~ 1 with a dispersion of approxi- 
mately 4 km s 1 . It is not clear what the origin of the discrepancy 
between these two cluster radial velocities is, but since the Jones 
et al. data are published, we adopt the median of their values 
as the cluster mean v r . Thus the radial velocity of HD 16605 is 
about 2cr from the cluster median. This discordant radial veloc- 
ity is possibly an indication of non-membership in NGC 1039, 
but with all other indicators of membership in favour of mem- 
bership, we still regard HD 16605 as a probable cluster member. 

5.2. HD 16728 = NGC 1039-307, B9V 

No previous magnetic observations of this star are available in 
the literature. We have two observations of the star. One has very 
low SNR for a reason that we have not been able to determine, 
and is omitted from Table [3] The second measurement yields 
a (B z ) value with an uncertainty of almost 300 G, mainly be- 
cause the LSD spectral line is shallow and broad, with a v e sin i 
of about 75 kms . The line profile shape could be that of an 
SB2. Synthesis of a small section of the ESPaDOnS / spectrum 
indicates that the atmospheric abundances in this star are close 
to solar values. (Systematic synthesis of ESPaDOnS / spectra 
will be a separate aspect of the investigation of these data which 
will be reported in a later paper. Here we simply provide a few 
first results that help us to separate magnetic Ap stars from other 
classes of stars.) It appears that this star is probably a normal star 
of r e ff ~ 1 1500 K rather than a magnetic Ap star. 

Proper motion data indicate that this star is probably a mem- 
ber of NGC 1039 (Baumgardt et al. [20001 Dias et al. |200TT l. 
Because it is not clear from our observations whether the star is 
single or is a spectroscopic binary, we cannot test cluster mem- 
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bership with a v r measurement. We consider HD 16728 to be a 
probable cluster member. 

5.3. HD 19805 = Melotte 20-167, B9.5V 

A large magnetic field of order 1 kG was reported by Bychkov et 
al. d2003 ), based on four unpublished measurements. However, 
three unpublished observations of the star using MuSiCoS, with 
ergs in the range of 125 - 175 G, do not reveal any sign of 
such a large field. Our two observations, with uncertainties for 
(B z ) of 24 and 14 G, also show no indication of any field, nor 
is a significant Zeeman effect detected in V. In addition, the 
spectrum of this star closely resembles that of a normal star of 
r eff * 10200 K, with strong He I and O I, normal Cr II, and 
no obvious sign of rare earths. Synthesis of / in small spectral 
regions confirms essentially solar abundances of these elements. 
We classify this star as normal. 

On the basis of proper motions, this star is a member of the 
a Per cluster, Melotte 20 (de Zeeuw et al. |1999l ). Robichon et al. 
( 1 19991 1 give the cluster velocity as -0.2 + 0.5 km s ; the values 
tabulated by Kharchenko et al. ( 12005b are similar. Our observed 
radial velocity v r = kms^'is consistent with membership; we 
consider this star to be a cluster member. 



5.4. HD 108945 = Melotte 111-160, A2p SrCr 

A weak field (109 ± 44 G) was definitely detected in this star by 
Shorlin et al. (2002). One measurement using FORS1 was re- 
ported in Kochukhov & Bagnulo (2006 ), where the field value is 
based only on the Balmer lines; no field was detected. In Table|4] 
we give the field value measured from these same ESO archival 
spectra using both Balmer and metal lines; again no field is de- 
tected, with cr B = 75 G. In two of our three ESPaDOnS measure- 
ments, (B z ) is non-zero at more than the 3cr# level, although the 
field is always of order 100 G or less. Furthermore, in all three 
measurements, a clearly non-zero signal is found in V, confirm- 
ing the reality of the small fields measured. HD 108945 is cer- 
tainly a magnetic Ap star. 

Proper motion data (cf Paper II) indicate that this star is a 
member of the Coma Berenices cluster. The mean cluster v r is 
-0.1 + 0.2 kms _I (Robichon et al. I1999t . Our measured radial 
velocity of +1 kms _I is consistent with membership; we con- 
sider this star a cluster member. 



5.5. HD 144661 = HIP 79031, B8p He-wk 

Three (B z ) observations of this He-weak star, with cr B ~ 200 G, 
were made by Borra et al. (Q983 ), but no strong evidence of a 
field was found. Our two magnetic measurements (<x B ~ 50 G) 
yield (B z ) values not significantly different from zero, and no 
sign of a magnetic signature is found in either V line profile. The 
LSD line profile is constant and appears to be a simple rotational 
profile. Borra et al. reported that this star is a member of the PGa 
class of He-weak stars, which seems to be a high-temperature ex- 
tension of the HgMn stars and, like the HgMn stars, seem to be 
chemically peculiar but non-magnetic. Woolf & Lambert ( 1999 ) 
suggest that the star is Hg-rich. Synthesis of a small section of 
spectrum shows that the atmosphere is depleted in He (by more 
than 1 dex) and S, almost normal in Si, and about 0.5 dex over- 
abundant in Fe (see also Norris [T971| ). We classify HD 144661 
as a non-magnetic, He-wk (probably PGa; cf Borra et al. 1983) 
star. 



On the basis of proper motions (de Zeeuw et al. 119991 ) the 
star is a member of Sco OB2. de Zeeuw et al. give a mean radial 
velocity of -4.6 km s for the Upper Sco region of the associa- 
tion. Our v r values are consistent with this value, and we consider 
the star an association member. 

5.6. HD 153948 = NGC 6281-15, Ap Si 

A field of about 200 G was detected in this star using FORS1 
(Paper I). All three of our measurements yield fields of about 
this size, although none of the (B z ) values differs from zero by 
3o"b. However, two of our observations yield definite detection 
of non-zero V signatures in the LSD lines. Furthermore, the line 
profile is strongly variable in shape, but not in v r . We classify 
HD 153948 as a definite magnetic Ap star. 

The proper motions are consistent with membership in 
NGC 628 1 (Paper II). The position in the HR Diagram is consis- 
tent with the cluster isochrone. On the other hand, the radial ve- 
locity measured on all three nights, about -2.5 kms -1 , is incon- 
sistent with the reported cluster radial velocity of -15.4 kms -1 
(Kharchenko et al. 2005 ). However, this value for the cluster ra- 
dial velocity appears to be based on only one measurement of 
one star, and is therefore highly uncertain. We disregard this in- 
consistency, and accept HD 153948 as a probable cluster mem- 
ber. 

5.7. HD 317857= NGC 6383-3, A1IVp 

Both FORS 1 (Paper I) and our single field measurement of this 
star reveal a strong (kG) field. HD 3 1 7857 is certainly a magnetic 
Ap star. However, both the radial velocity of -24 kms -1 mea- 
sured (once) by Lloyd Evans ( 19781 ), and the value of - 10 km s _1 
from our ESPaDOnS spectrum, are inconsistent with the clus- 
ter radial velocity of about +5 km s reported by Kharchenko et 
al. (|2005 ); although the two stellar measurements taken together 
suggest that the star is a spectroscopic binary. More definitively, 
HD 317857 is about 0.8 dex too bright for the cluster isochrone. 
It appears that this star is a probable cluster non-member. 

5.8. HD 318100= NGC 6405-19, B9p 

Both the FORS 1 measurement and the present ESPaDOnS ob- 
servation reveal longitudinal fields of a few hundred G, signifi- 
cant at the several <x B level. The star also has an extremely large 
value of the Aa peculiarity parameter (Maitzen & Schneider 
rT984b . HD 318100 is definitely a magnetic Ap. 

Kharchenko et al. (2005) report proper motions that are 
consistent with the cluster mean values; our radial velocity 
of -8 kms -1 is consistent with the cluster mean of -7 + 
1 km s (Kharchenko et al. 120051 ); and the inferred luminosity 
(Paper II) is consistent with the cluster isochrone. We consider 
HD 3 1 8 1 00 a member of NGC 6405 . 

5.9. HD 320764 = NGC 6475-23, A1Van 

One FORS1 measurement, based mainly on Balmer lines, 
yielded no field detection with erg = 66 G. Our two measure- 
ments have uncertainties of the order of 4-500 G due to the ex- 
tremely broad and shallow LSD line, and provide no really use- 
ful constraint. Folsom et al. (2007 ) report v e sin i - 225 km s , 
and their abundance analysis shows that this is a normal star, not 
a magnetic Ap. 
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Proper motions suggest that this star is a probable cluster 
member (cf Paper II). 

5.10. HD 162576 = NGC 6475-55, B9.5p SiCr 

Our three field measurements with ESPaDOnS are the first made 
of this star. In spite of standard errors as low as about 20 G, none 
of the measurements of (B z ) are significantly different from zero. 
However, we have one significant detection of a non-zero V sig- 
nature. In addition, the LSD profile is strongly variable, without 
any change in v r . The star has been reported to be photometri- 
cally variable (North 1 1984b . and Folsom et al. (120071) report that 
the abundances found in the atmosphere are clearly typical of 
magnetic Ap stars. We classify HD 162576 as a magnetic Ap 
star; it appears to be close to the weak-field limit of such stars 
(cf Auriereet al. 120071 

Astrometric data indicate that the star is a member of 
NGC 6475 (cf Paper II), and our measured values of v r are con- 
sistent with the cluster radial velocity of -14.7 + 0.2 kms"' 
(Robichon et al. 1999). This star is a probable cluster member. 

5.11. HD 162588 = NGC 6475-59, Ap Cr 

None of our three field measurements (the first such observations 
of this star) reveal a (B z ) value that is significantly different from 
zero, but the <x B s are fairly large, of order 100 G, because of the 
elevated value of v e sin i (70 km s~ ). No significant V signal is 
detected. However, the LSD / line profile is strongly variable 
without significant change in v r . The star is a photometric vari- 
able (North n~9 841 ). We have synthesised a small section of the 
/ spectrum, and find that Cr is overabundant by about 1.5 dex. 
In spite of the absence of direct evidence for the presence of a 
magnetic field, we classify the star as a magnetic Ap star. 

Astrometric data indicate that HD 162588 is a member of 
NGC 6475 (Dias et al. 2001). Our measured v r values are con- 
sistent with the cluster velocity, and we consider the star a cluster 
member. 

5.12. HD 162630 = NGC 6475-63, B9V 

No significant (B z ) field was detected in two measurements, and 
no significant V signature is present. The LSD profile is consid- 
erably deeper using a solar line list than with an Ap line list of the 
same temperature, which is usually a symptom that many lines 
present in the Ap list but not in the solar list are actually weak 
or absent in the analysed star. In addition, the LSD line profile is 
apparently a simple rotation profile, and it does not vary between 
our two observations. In the / spectrum, He I and O I are quite 
strong, which is normally not the case for magnetic Ap stars. 
Synthesis of a small section of / spectrum shows that Fe appears 
underabundant relative to solar by about 0.3 dex, while Si and He 
have abundances close to solar. We conclude that HD 162630 is 
not a magnetic Ap star, but an (approximately?) normal star of 
r eff ~ 10400K. 

HD 162630 does appear, both on the basis of proper motions 
and of our v r measurements, to be a member of NGC 6475. 

5.13. HD 162656 = NGC 6475-72, B9V 

Our single measurement, the first field measurement of 
HD 162656, shows no sign of a magnetic field either in the value 
of (B z ) or in the LSD V profile. The LSD profile is deeper, and 
the field measurement more accurate, using a solar line list than 



with an Ap list. The LSD line profile appears to be broadened 
simply by rotation. Although the star has a T e ff value of about 
9500 K, the / spectrum shows strong He I and O I lines, and the 
Cr II lines near 6150 A are weak or absent. Synthesis of short 
sections of the / spectrum confirms these qualitative observa- 
tions; He, Si, and Fe all have abundances close to solar, while Cr 
appears to be about 0.2 dex below solar. We classify this star as 
normal. 

The proper motions are consistent with membership in 
NGC 6475, but our one measured v r value is not. However, it 
has been reported by Gieseking (119771 ) that the star is a spectro- 
scopic binary, and that the y value of the orbit is consistent with 
cluster membership. We consider the star to be a cluster member. 



5. 14. HD 162725 = NGC 6475-88, AOp SiCr 

Our data are the first magnetic field measurements of this star. 
One of our two measurements reveals a (B z ) value lcr B different 
from zero, although the other is consistent with zero. Both ob- 
servations have LSD V signatures that are clearly different from 
zero. Folsom et al. (2007 ) report that the atmospheric composi- 
tion of HD 162725 is clearly that of an Ap star. The star is also 
a photometric variable (North [T984I ). It is clear that this star is a 
magnetic Ap star. 

Both proper motions and our radial velocity are consistent 
with classification of the star as as member of NGC 6475. 



5. 15. HD 169842 = NGC 6633-39, Alp SrCr 

Kudryavtsev et al. (2006) report probable detection of a field 
of B rms « 370 G from five measurements with typical standard 
errors of about 180 G. Our three observations provide clear de- 
tections of (B z ) at about the l0cr B level, and confirm the general 
magnitude of the field reported by Kudryavtsev et al. Synthesis 
of two spectral windows (cf. Figure [TJ indicates that Cr is over- 
abundant relative to solar abundances by about 2 dex, while Si 
and Fe are near solar abundances. HD 169842 is definitely a 
magnetic Ap star. 

Proper motions are consistent with membership in 
NGC 6633. The cluster radial velocity is reported by 
Kharchenko et al. (120051 ) to be between -25 and -29 kms . 
Individual radial velocities from references cited in WEBDA 
support a value around -29 km s~ . Our measured v, values are 
consistent with these values, and we consider the star to be a 
cluster member. 



5.16. HD 169959 A = NGC 6633-58, AOp Si 

The field of this star was first detected with FORS1 (Paper I). 
Our field measurement confirms the presence of a field and its 
general magnitude in this magnetic Ap star. 

The proper motions and previous radial velocity measure- 
ments (Hintzen et al. I1974I ) are consistent with membership in 
NGC 6633. However, our one radial velocity (-53 km s~') is in- 
consistent with both the cluster velocity of about -30 km/s and 
with the measurements of Hintzen et al. It is not at all clear what 
the origin of this discrepancy is. Furthermore, HD 169959A, if it 
is a cluster member, is a (very) blue straggler. We conclude that 
it is a probable cluster non-member. 
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5. 17. HD 1 70054 = NGC 6633-77, B6IV 

No field was detected in two previous FORS1 measurements 
with standard errors of about 75 G (Paper I). Our two new mea- 
surements, one with a larger <jb (due to a period of low trans- 
parency) and one with smaller cr B also reveal no field, nor is 
there any hint of a field in the V signatures. Spectral lines of 
He I, O I, and Ne I are strong, while Cr and Fe appear roughly 
normal. Synthesis confirms abundances of He and Fe that are 
close to solar, while Si appears about 0.2 dex below solar. We 
classify HD 170054 as a normal B star, although the peculiar 
shapes of the / line profiles, and their variability, suggest that 
this star may be a pulsating variable. 

Proper motions indicate that the star is a cluster member, and 
our v r value is consistent with this. We consider the star to be a 
cluster member. 

5. 18. HD 1 70860 = IC 4725-153, B9IV/Vp 

The magnetic field of this star was not detected in one 
FORS1 measurement, which achieved <x B = 69 G. Our single 
ESPaDOnS measurement nominally provides about a 2<tb de- 
tection, although both (B z ) and <x B are not very well determined, 
due to uncertainties as to where the edges of the LSD line pro- 
file are. However, a large and clearly significant V signature is 
present. In addition, the Aa value (Maitzen 1985) is rather large, 
providing further support to our classification of this star as a 
magnetic Ap star. 

Proper motions indicate that HD 170860 is a mem- 
ber of IC 4725. The mean cluster radial velocity has been 
reported as having values from -4 + 4 (Feast 1 19571 1 to 
+2.4 kms _1 (Kharchenko et al. 12005b . Because of the difficulty 
of determining the centroid of the LSD line, our v r value is not 
derived from this line; instead, it is derived by fitting Gaussians 
to the cores of the Balmer lines Ha - H<5, and the strong 
Si II lines 4128-30 and 6347-71 A. The value found by us, 
-13 ± 3 kms -1 , is not the same as the cluster mean, but with 
only one measurement, we cannot rule out the possibility that 
the star is a spectroscopic binary. We regard HD 170860 as a 
probable cluster member. 

5.19. HD 1 72271 = IC 4756-1 18, AOp Cr 

Our measurement is the first reported magnetic observation of 
this star. We find a marginally significant (B z ) in this star. The 
significantly peculiar value of the Geneva index Z = -0.027, 
and the apparent absence of He I lines in HD 172271, strongly 
suggest classification as a magnetic Ap star. This is confirmed 
by synthesis of a small spectral region, from which we find that 
Cr is overabundant by about 2.5 dex, while Fe is about 0.5 dex 
overabundant. 

The value of fi a cos 6 is about 2<x from the cluster mean, 
but yU (5 is consistent with membership. The cluster v, is -25.8 + 
0.2 kms 1 (Robichon et al. 1999; note sign error on their v r 
value); our measured v r = -21 ±2-3 kms _1 (the star has a 
large v e sin z) is reasonably consistent with this value. We regard 
this star as a probable cluster member. 

5.20. HD 205073 = NGC 7092-69, Al 

Our two measurements of this star are the first reported mag- 
netic observations. No field is detected, either in the (B z ) values 
or in the V profiles. Considering the rather small v e sin i value 
(16 kms _I ), the lack of detection is strong evidence that this 



star is not a magnetic Ap. Furthermore, the field measurement is 
more precise using a solar composition line list for LSD analy- 
sis. Neither the Aa nor the Geneva Z peculiarity indices are sig- 
nificantly different from zero, and our spectra show line profiles 
which appear to be simply rotationally broadened, with clearly 
present He I and O I but very weak Cr II and no rare earths. 
The abundances of Cr and Fe deduced from spectrum synthesis 
of small sections of the spectrum are close to solar. It is clear 
that this star is a normal A star. Note incidentally the two quite 
different v r values of the two observations; HD 205073 is a spec- 
troscopic binary. Since close binaries are rare among magnetic 
Ap stars, the radial velocity variability supports our view that 
HD 205073 is not such a star. 

The proper motions are consistent with cluster membership 
(Baumgardt et al. 2000; Dias et al. 1200 11 1. We consider the star 
to be a cluster member. 



5.21. HD 205331 = NGC 7092-118, A1V 

The first published magnetic measurements of this star were 
those of Kudryavtsev et al. {2006 ), who failed to detect a field 
in three observations with cr B of about 280 G. Our three data re- 
duce the standard error to the 30 - 40 G range. We do not have 
a statistically significant detection in the V signature in any of 
the three observations, but one (B z ) value (integrated over the 
whole broad line) is about 3<tb different from zero, and suggests 
a field of about 140 G. However, our field measurement is as 
precise with a solar line list as with an Ap list, and the LSD pro- 
file appears to be due to pure rotational broadening. Both these 
features are often signatures of a normal star. Our spectra show a 
clear He I /15875 line (in spite of the low T e s value), strong O I, 
weak Cr II, and no obvious rare earths. The brightness is con- 
stant in the Hipparcos photometry catalogue. Synthesis of two 
short spectral regions in the ESPaDOnS / spectrum yields near- 
solar abundances of Si, Cr and Fe, assuming that the microtur- 
bulence parameter is about 2 km s . We classify HD 20533 1 as 
a normal star of T e s ss 9450 K in spite of the single apparently 
significant field detection. 

The proper motions and parallax are consistent with cluster 
membership (Baumgardt et al. 120001 Dias et al. |2001| l. The ac- 
curate cluster radial velocity of -5.4 ± 0.4 km s (Robichon et 
al. [1999) is agrees with the three values of v r measured by us. 
We consider the star a cluster member. 



5.22. BD +49 3789 = NGC 7243-490, B7p Si 

Our single measurement of the field of this star is the only one 
available. The precision is very low (cr B ~ 400 G), and no detec- 
tion is achieved. The large values of the photometric peculiarity 
indicators Aa and Geneva Z, and the presence of strong Si II 
lines in our / spectrum suggest that BD +49 3796 is actually a 
magnetic Ap star. This view is confirmed by synthesis of two 
small spectral windows, which yields an underabundance of He 
by about 1 dex, and overabundances of Si, Cr and Fe by respec- 
tively 1, 0.7, and 0.5 dex relative to solar abundances. 

Astrometric data for this star, and HR Diagram position, are 
consistent with cluster membership, but our single v, value dif- 
fers from the cluster mean of between -9 km s -1 (Hill & Barnes 
[T97fl > and -13 kms" 1 (Kharchenko et al. [20051 ) by about 5 
km s _1 . However, Hill & Barnes report four v r measurements be- 
tween -13 and +32 km s , and consider the star a spectroscopic 
binary. We regard the star as a probable member. 
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5.23. HIP 10991 1 = NGC 7243-370, AOVp Si 

Our magnetic measurement is the only one available for this star, 
and does not reveal a significant (B z ) or V signal, although there 
is a hint in our low-precision 2<x B (B : ) that there might be a field 
of a few hundred G present. In any case, the photometric pecu- 
liarity indices Aa and Geneva Z are large enough to strongly sug- 
gest that HIP 10991 1 is a magnetic Ap, and the obvious strength 
of Si II and Cr II lines supports this view. Synthesis of a small 
spectral window yields He depleted by 1 .5 dex, and Si and Fe en- 
hanced by 1 and 1 .2 dex respectively relative to solar. We classify 
the star as a magnetic Ap star. 

The proper motions (Baumgardt et al. 2000; Dias et al. 2001 ) 
and position in the HR Diagram are consistent with cluster mem- 
bership. Hill & Barnes ( 119711 ) report four v r values between -8 
and -24 kms -1 but do not flag the star as a velocity variable, 
perhaps reflecting the relatively large mean errors. Our single v r 
value, with an uncertainty of roughly 2 km s , is consistent with 
membership. We consider the star a cluster member. 

6. Characteristics of the data set 

In three nights (of which all but a few hours were clear) with 
ESPaDOnS we have obtained 44 observations of 23 stars which 
are possible magnetic Ap stars and possible members of open 
clusters. For 12 of these 23 stars, our measurements represent 
the first published magnetic data. 

Because of the time of year when our observing run was 
scheduled, about half of our ESPaDOnS observations were made 
on stars which had been previously observed, either earlier in our 
cluster survey with FORS 1 (Paper I), or by others. Nevertheless, 
the data we have acquired provide valuable new information 
about possible cluster magnetic stars, even when these stars have 
been previously measured for magnetic fields. 

One important class of information provided by our new 
ESPaDOnS data is accurate radial velocities of almost all the 
stars observed. Because the usefulness of proper motions as a 
discriminant of cluster membership diminishes roughly propor- 
tionally to distance (the required accuracy of motions increases 
with distance, but the uncertainties of available data set remain 
roughly constant), while radial velocities continue to be a pre- 
cise discriminant to whatever distance they can be measured ac- 
curately, the radial velocities are particularly useful in testing 
membership of stars in clusters that are more than 1 - 200 pc dis- 
tant. For all but a few of the 23 stars observed, our data provide 
a valuable test of membership, mostly supporting cluster mem- 
bership. However, our v r data indicate that two stars previously 
thought to be cluster members (HD 317857 and HD 169959 A) 
are probable non-members. 

Another important type of information provided by 
ESPaDOnS is clearer discrimination between stars which are 
magnetic Ap stars with fields small enough to be difficult to de- 
tect, and stars which are not magnetic Ap stars. In order to obtain 
representative field strength distributions for various star masses 
and ages, we have included a number of stars in which fields 
have not yet been detected, but for which other evidence exists 
indicating that they are probably magnetic Ap stars, in the sam- 
ple analysed in Paper II. Even with ESPaDOnS observations of 
such stars it is not always possible to detect or rule out a field, 
but in most cases the ESPaDOnS spectrum allows one to classify 
a star with confidence as either a magnetic Ap or a normal star, 
something which is not usually possible with FORS 1 / spectra 
(see for example the discussions of HD 162630, HD 162656, 
HD 170054, HD 205073, and HD 205331). Our / and V spec- 



tra allow us to categorise nine stars of the present sample as not 
magnetic Ap stars; in most cases, these appear to be normal stars. 

One obvious aspect of the present data set is that in stars 
with relatively rich and sharp-lined metal spectra, ESPaDOnS 
data may allow the detection of fields not readily detected with 
FORS 1 . Of the two stars for which we report the first direct evi- 
dence of a magnetic field, one (HD 162725) was previously ob- 
served with FORS 1 (Paper I). 

Of course, for all the stars we observe, the new measure- 
ments contribute to convergence of B rms , the RMS longitudinal 
field (see Paper II and below), which is the measure we use of 
the global characteristic field strength of each star in our sample. 



7. Magnetic field evolution through the main 
sequence 

We next use the new data acquired from ESPaDOnS to revise 
our global database of magnetic Ap stars that are members of 
clusters. We flag stars that are found to be probable or certain 
non-members, or probable or certain non-magnetic stars. We add 
newly studied probable or certain magnetic cluster Ap stars to 
the database. For each of the new stars, we determine the funda- 
mental parameters T e g and log(L/L ) using the same techniques 
as in Paper II. As discussed in detail in that paper, we assume 
that the uncertainty in log(r er r) is ±0.02, and the uncertainty in 
log(L/L ) is +0.1. With these data, we place each star in the 
cluster HR diagram, check whether its position in the HR dia- 
gram is close enough to the cluster isochrone to be consistent 
with membership, and estimate its mass from its position. 

For each retained magnetic cluster member (old or new) in 
the database for which we have new magnetic data, we compute 
the revised characteristic field strength B rms , which is simply the 
RMS average of all available field measurements, or the subset 
with the smallest uncertainties. (Although B rms is a fairly crude 
measure of stellar field strength, it is the best measure we can 
easily form from the data until enough is known about most stars 
in the sample to actually compute simple models of individual 
field structure.) Finally, we re-examine the evolution of magnetic 
fields with time elapsed on the main sequence to see if any of the 
conclusions of Paper II have to be modified. 

After carrying out these steps, we find that we have 12 stars 
which satisfy our criteria of cluster membership and probable or 
definite magnetic Ap nature, and for which the magnetic field 
data are new or improved, so that a new value of B ims may be 
calculated. These stars are listed in Table [5] The table is similar 
to Table 3 of Paper II in the information it contains. For each 
magnetic Ap cluster member for which new magnetic data have 
been obtained, the table lists the cluster to which the star be- 
longs, the cluster age and uncertainty (see below), the true (not 
apparent) distance modulus DM of the cluster, the star identi- 
fication, the star's log(r e ff), log(L/L ), M/M Q with uncertainty, 
fractional age t with uncertainty, the revised value of B rms , and 
the number ng of magnetic observations contributing to B rms . 
(Note that for Ap stars for which we have not yet detected a 
field, we use only the most accurate measurements in forming 
firms, so as not to inflate that parameter by including numbers 
with unnecessarily large uncertainties.) 

We also use the newly measured stars to test the cluster ages 
found in the literature by comparing the position of stars near 
the terminal age main sequence (TAMS) with the isochrones. 
In Paper II, we found that this procedure made it possible to im- 
prove the precision of age determination of several open clusters. 
Compared to Paper II, we obtain improved ages with smaller er- 
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ror bars for NGC 6475, IC 4725, and IC 4756. These improved 
ages are listed in column 2 of Table [5] 

Then finally we include the new data together with the data 
from Paper II to check whether the new magnetic measurements 
alter our previous conclusions (Paper II) concerning the evo- 
lution of magnetic fields with time through the main sequence 
lifetime of magnetic Ap stars. The full, revised data set of stars 
which are probable cluster members, probable Ap stars, and for 
which at least one field measurement is available, is used to cre- 
ate FigureHJ which contains information similar to that found in 
Figure 5 of Paper II. 

In this figure the RMS field estimates fi rms for individual 
stars (in three mass bins: 2-3 M/M e , 3-4 M/M Q , and 4 - 
5 M/Mq, bins for which a significant number of stars are now 
available) are plotted as a function of cluster (and thus of stellar) 
age (left three panels) and of fractional age (right three panels). 
The full sample is divided into mass bins in order to be able to 
study the variation of typical field strength fi rms with stellar age 
as a function of stellar mass. (Note that the star NGC 2244-334, 
with firms = 9515 G, log(age) = 6.40, and t = 0.02 is off-scale 
in the upper left corner of both bottom panels.) 

The present sample differs from the sample of Paper II essen- 
tially in that the first magnetic measurements for five stars have 
been added, and the fi rms data for five other stars have changed, 
always by a factor of less than two. Although the ESPaDOnS 
data do not increase the total size of the sample very much, these 
observations have substantially increased the information avail- 
able about lower-mass Ap stars that have completed a large part 
of their main sequence evolution (see the fractional age column 
of Table 01. These data provide a significant improvement for a 
region of parameter space that was not very well represented in 
the previous sample. 

The data of the full sample of probable cluster magnetic Ap 
stars are displayed here a little differently than in Paper II. In 
Figure|4]we show the variation of fi lms (rather than its logarithm) 
as a function of stellar (= cluster) age in the left panels, and as 
a function of fractional age in the right panels. The three panels 
in each column show the data for the three mass bins mentioned 
above. In all three right-hand panels we can now see clearly that 
RMS fields larger than 1 kG are common at small fractional age, 
while the (now numerous) stars of fractional ages larger than 
about 0.4 for 2 - 3 M , 0.2 for 3 - 4 M , and 0. 1 for 4 - 5 M 
have fields that are always less than 1 kG. These trends are not 
so easily seen in the left panels using log(age) as the abscissa, 
which expand the region of small fractional age at the expense of 
regions of large fractional age, but close examination of the left 
panels shows that the same effect is present. In all three panels 
on the left side of the figure, a significant number of low field 
stars (but no high-field stars) are present near the TAMS. 

The obvious interpretation of this figure is that the RMS field 
declines rather sharply with age after a characteristic time which 
is of the order of 250 Myr for stars in the mass range 2-3 M , 40 
Myr for 3 - 4 M , and 15 Myr for 4 - 5 M stars. It is remarkable 
that the characteristic time of the field decreases rapidly as the 
stellar mass increases, both in absolute value and as a fraction of 
the stellar main sequence lifetime. 

One expected physical effect which would lead to some sub- 
stantial decline in surface magnetic field is of course the geo- 
metric expansion of the stars by a factor of about two in radius 
as they evolve from the ZAMS to the TAMS. If magnetic flux 
(and global field topology) are approximately conserved, the in- 
crease of radius by a factor of order two will lead to decrease in 
field strength by a factor of order four. We test this possibility by 
plotting fi rms x (R/R Q ) 2 for the stars of Figure|4]as a function of 



log(age) and t in Figure [5] This quantity is a reasonable proxy 
for the emergent or apparent magnetic flux (normalised to an av- 
erage field of 1 G over a star of one R Q ); it differs from the true 
normalised emergent flux by a factor (dependent on the view- 
ing geometry and obliquity of the stellar field) which relates the 
RMS value of the observed sample of values of mean longitudi- 
nal field (B z ) to the typical field modulus at the stellar surface. 
This factor is at least about 3, and is likely to be of the order of 5 
for most of the stars of our sample. Of course, the emergent flux 
may or may not be closely related to the flux deep in the stellar 
interior. However, note that fields are detected among the older 
stars in all our mass bins. The fields decline in all bins, but not 
to zero in the main sequence lifetime of the stars considered. 

Overall we see that Figure|5]shows much the same behaviour 
as Figure [4] Near r = there is a rather large range of es- 
timated normalised emergent magnetic fluxes, and a few stars 
which exceed the median normalised flux value by a factor of 10 
or more. As one moves towards larger fractional age, the fluxes 
converge on lower values - quickly for the higher mass ranges, 
more slowly for the (top) 2-3 M bin. From the right column of 
figures, one plausible hypothesis is that stars with initially low 
fields, less than about 1 kG, tend to conserve flux as they evolve 
across the main sequence. However, the stars with initially large 
fields, those that contribute the large normalised flux values in 
the figure, appear to lose much of this flux within a time short 
compared to the main sequence lifetime. This flux loss is not 
due to the geometric expansion of the star, already taken into ac- 
count in the method of estimating the normalised magnetic flux. 
Instead, it could be some kind of relaxation or redistribution pro- 
cess 

We do not find any analogues of the young, high-field (and 
high emergent flux) stars among the (now reasonably numerous) 
older stars in any of the mass bins considered here. We conclude 
that the present data may be consistent with approximate emer- 
gent flux conservation in most magnetic Ap stars, but that the 
initially high-field stars seem to suffer a real decline in the to- 
tal emergent magnetic flux during the main sequence phase of 
evolution. Alternatively, the present data may show that typical 
emergent magnetic flux declines significantly in most or all stars 
(rather than only in the stars with initially large fields and fluxes), 
by a factor of two or three, during the main sequence phase. 

It appears that our observations are consistent with the fossil 
field hypothesis (as we have implicitly assumed in discussing the 
possibility of flux conservation above), in that the flux in older 
stars are within a factor of two or three of the fluxes in young 
stars in each mass bin. As many authors have noted, these fluxes 
do not show the expected correlation with rotation that is found 
in the presumably dynamo-generated fields of lower main se- 
quence stars. In fact, some of the largest fields are found among 
stars that rotate very slowly. 

Furthermore, one theoretical concern about the fossil field 
hypothesis has be the question of whether fossil fields are stable 
over long time scales. Tayler (U973I I, Markey & Tayler ( 119731 
1 1974| l, and Wright (U973I I have shown that fields which are 
purely poloidal or purely toroidal in the radiative interior of a 
star are unstable on time-scales that are short compared to the 
main sequence lifetime of an Ap star. The fact that we find fields 
with similar fluxes at all ages in each mass bin certainly suggests 
that, if these fields are fossils, nature has been able to resolve this 
problem. 

This problem has recently been clarified on the basis of a 
set of important calculations discussed by Braithwaite & Spruit 
(120041) and Braithwaite & Nordlund ( 120061) . This study has been 
primarily focussed on the issue of the long-term stability of fos- 
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Fig. 4. This figure shows the currently estimated values of B rms as functions of logarithmic stellar age (left) and of fractional age 
(right) for three mass bins, from top to bottom 2 - 3M , 3 - 4M , and 4 - 5M . Filled symbols are stars for which a field is definitely 
detected; open symbols are probable magnetic Ap stars in which no field has yet been detected. The right-hand limit of each of the 
panels using log(age) as abscissa is near the main sequence lifetime for stars in that mass range. In the bottom pair of panels, one 
point (for NGC 2244-334) has such a large field (B rms = 9.52 kG) that it is off scale (at log t = 6.4 + 0.10 and r = 0.02 + 0.01 
respectively). 



sil magnetic fields in Ap stars (and magnetic white dwarfs). 
These authors conclude, on the basis of numerical simulations 
of field evolution in a non-rotating polytropic star, that the insta- 
bilities found in earlier work lead to rapid field evolution which 
generally creates a mixed poloidal-toroidal field structure with 
a twisted toroidal field. Once formed, this structure appears to 
be stable on a longer time-scale. This work appears to provide 
confirmation that stable fossil field structures exist, and can be 
achieved by a star starting from a variety of initial fields, in 
agreement with the observations of fields with comparable fluxes 
(to within plus or minus an order of magnitude) in Ap stars of 
all main sequence ages. 

Since this work studies the evolution of field structure in ra- 
diative stars similar to main sequence stars, it also leads to pre- 
dictions of how a global fossil magnetic field should evolve with 
time, at least under the hypotheses adopted in the model. The 
evolution found by Braithwaite and collaborators, when their 
computations are rescaled to approximately realistic magnetic 
diffusivities, is that the emergent flux gradually grows with time, 
by a factor that may be of order ten, until finally the main toroidal 
flux tube begins to emerge from the star, at which point the sur- 
face field declines rapidly. (The emergent flux is mostly poloidal, 
since a stellar atmosphere does not support the large current den- 
sities needed for a strong surficial toroidal field component.) 
This evolution is estimated (roughly) to take somewhat longer 
than the main sequence lifetime of most Ap stars. Thus these 



computations predict that we should mainly see stars in the phase 
of surface field strength increasing with time. 

This is, of course, the opposite of the situation revealed 
by our data, in which we see both field and flux decline with 
time, at least in some stars and perhaps in all. Thus it appears 
that the field evolution predicted by the numerical modelling of 
Braithwaite and collaborators is not fully consistent with obser- 
vations. The observations are consistent in that fields appear to 
be stable on a long enough time scale to persist throughout the 
main sequence phase, in agreement with the numerical models, 
but the observed and computed fields appear to have somewhat 
different time evolution. 

However, it is perhaps not very surprising that the field evo- 
lution predicted by these numerical models differs from what is 
observed. This modelling neglects two effects that may signifi- 
cantly affect the evolution of surface fields. First, the computa- 
tions are for non-rotating stars. The effect of reasonably rapid 
rotation (and most magnetic Aps, including most of the stars in 
our study, have rotational periods of a few days, and equatorial 
velocities of several tens of km s~ 1 ) is to produce meridional cir- 
culation. In turn, this circulation transports angular momentum 
and almost certainly produces shear within the star. This shear 
will certainly strongly distort the internal field, which will prob- 
ably substantially alter the emergent flux. Secondly, the mod- 
elling has been carried out for a star of constant size and struc- 
ture. During the main sequence phase, a real star changes its 
overall structure significantly, and in particular increases in ra- 
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as Figure [4] 



dius by a factor of order two. This radius increase will counteract 
to some extent the increases surface field strength (although not 
the increased surface flux). It appears that further computations 
of field evolution, with incorporation of additional physics, will 
be of great interest. 

8. Conclusions 

In order to study the evolution of magnetic field strength and of 
atmospheric chemistry in magnetic Ap stars through their main 
sequence lifetimes, we have been carrying out a large survey 
of such stars in open clusters. As discussed in Papers I and II, 
the main interest of such a sample is the fact that the ages of 
such stars, especially during the first half of their main sequence 
lifetime, can be determined with considerably greater precision 
than is the case for magnetic Ap stars in the nearby field. 

Most of the survey so far has been carried out using FORS 1 
at the ESO VLT (cf Paper I), with an important contribution from 
earlier studies of the Ori OBI and Sco OB2 associations (see 
for example Borra [1981] and Thompson et al. 119871 ). In this pa- 
per we discuss the first use of the new instrument ESPaDOnS at 
the Canada-France-Hawaii Telescope for this survey. The three 
nights of data so far obtained from this instrument have made 
an important additional contribution to the overall survey. They 
have also made possible several very important checks of previ- 
ous data, and provided kinds of information not readily available 
from the lower resolution FORS 1 spectra. 

Comparison of ESPaDOnS and FORS1 measurements of 
the same stars show that for stars that are not too hot (say be- 
low about 12-13 000 K) and that have v e sin i less than roughly 



60 kms , the precision in (B z ) obtainable with ESPaDOnS in 
a given time on target is comparable to, and sometimes sub- 
stantially better than, the precision obtained with FORS1, even 
though the CFHT has only about 20% of the collection area 
of one VLT unit telescope. This is possible essentially because 
ESPaDOnS has far higher resolving power (65 000 compared to 
2 000 or less) than FORS 1 , and thus can fully exploit the Zeeman 
signal contained in the metal lines. If the star is rotating slowly 
enough, and is cool enough, this gives ESPaDOnS a very large 
advantage over FORS 1 observations, which must rely to a large 
extent on the signal in the Balmer lines. 

In addition, with ESPaDOnS we have the possibility of de- 
tecting a real and significant circular polarisation V signal in the 
averaged line profile even if the mean longitudinal field (B z ) is 
close to zero. In weak-field stars such as HD 108945 the field is 
much easier to detect this way than through conventional mea- 
surements of (B z ). This capability is also very helpful in decid- 
ing which proposed Ap stars in a cluster are in fact Ap stars, and 
which are normal or belong to other, non-magnetic, peculiarity 
classes. Having clear evidence on this point is very helpful in 
cleaning the field measurement sample of non-magnetic stars. 

Furthermore, with ESPaDOnS we obtain not only spectropo- 
larimetry, but also high-resolution / spectra over a wide wave- 
length window. These data are very useful for deciding which 
stars have the chemical peculiarities of magnetic Ap stars, and 
which appear chemically more or less normal. We have exploited 
this possibility above to eliminate some stars that have been sug- 
gested as magnetic Ap stars from the survey sample. In a later 
stage of this project, we plan to model in detail a number of the 
/ spectra from ESPaDOnS with the goal of obtaining a clearer 
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picture of the evolution of atmospheric chemical peculiarities 
during the main sequence lifetime of such stars. This cannot be 
done at present with useful accuracy for the low resolution spec- 
tra obtained with FORS 1 . 

The accurate radial velocities that can be obtained with 
ESPaDOnS (recall that the precision of radial velocities, for stars 
with v e sim of the order of 50 kms -1 , sometimes with non- 
symmetrical LSD profiles, is roughly ±1-2 kms 1 ) provide 
valuable tests of cluster membership for stars of our sample. The 
accuracy of radial velocity measurements from our FORS 1 data 
is probably not good enough to provide really useful constraints. 

The precision of the best field measurements obtained with 
ESPaDOnS is extremely high. The best measurement have stan- 
dard errors in the range of 15-30 G, as small as or smaller than 
the highest precision (about 30 G) that has been obtained reliably 
with FORS 1 . This has made it possible to add five new stars to 
the sample of cluster Ap stars with field measurements, and to 
add further measurements of several others stars of the sample. 
The new measurements provide valuable additional data on stars 
having low mass and large fractional ages t, a part of the sample 
which was under-represented in our previous data set. 

We have re-examined the full sample, and present the data in 
a somewhat different form than in Paper II. We find that our ear- 
lier conclusions (see Paper II) are somewhat altered. As before, 
we find that magnetic fields of stars with masses above about 
3M decline rather strongly through the main sequence life, 
dropping in strength by a factor of several after about 4 10 7 yr 
and 1.5 10 7 yr for stars in the 3-4 and 4-5 M bins respec- 
tively. This is consistent with the 3 10 7 yr found for stars above 
3M in Paper II. However, in contrast to the results of Paper II, 
we now find clear evidence that the fields of lower mass Ap stars 
(2-3 M ) appear to decline with increasing age after roughly 
2.5 10 8 yr. Furthermore, for all masses, the data suggest that even 
the emergent magnetic^Mx either declines somewhat in all stars 
through the main sequence lifetime, or that the emergent flux in 
the more strongly magnetic stars found at young fractional ages 
is somehow reduced. 

When our results are compared with the numerical computa- 
tions of Brathwaite & Spruit ( 2004 ) and Braithwaite & Nordlund 
(2006 ), the data are found to support the numerical results in that 
both observations and numerical simulations agree that the ap- 
parently fossil fields of magnetic Ap stars are stable on an evo- 
lutionary time scale. On the other hand, the numerical models 
predict that during the main sequence phase the typical mag- 
netic fields of Ap stars should be observed to increase, the op- 
posite of what is observed. However, since these computations 
neglect both the internal (poloidal and toroidal) circulation flows 
induced by rotation (which will probably strongly distort the 
field), and the global radius increase as the star evolves from 
ZAMS to TAMS, it is quite possible that the evolutionary predic- 
tions will be substantially modified as more physics is included 
in the modelling. It will be of great interest to see whether in- 
cluding these (or other) physical effects in the numerical models 
leads to more complete agreement between computation and ob- 
servation. 

We are continuing our survey, both to obtain more precise 
values of the characteristic field strength indicator B lms used in 
this study for stars already observed, and to increase the sam- 
ple to a larger number of clusters and stars. Our ultimate goal is 
to provide clear and unambiguous observational constraints on 
the evolution of both magnetic field strength (and perhaps even 
structure), and of atmospheric chemistry, through the main se- 
quence lifetime of the magnetic Ap and Bp stars. 
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Table 3. Magnetic field measurements of possible cluster Ap stars 



Star 


HJD 


SNR 


v r 

(kms" 1 ) 


<5z> ± 

(G) 


fid 


HD 16605 


2453570.110 


608 


-1 


-2099 ± 34 


DD 




2453571.081 


689 


-1 


-2205 ± 29 


DD 


HD 16728 


2453572.123 


1510 




-2 ± 272 




HD 19805 


2453570.131 


1045 





6 ±24 






2453572.129 


1453 





-3 ± 14 




HD 108945 


2453569.732 


1884 


+1 


-111 ±28 


DD 




2453569.740 


1872 


+1 


-75 ± 29 


DD 




2453570.735 


2487 


+1 


109 ± 20 


DD 


HD 144661 


2453571.758 


2814 


-3 


-55 ± 48 






2453571.758 


2946 


-4 


11 ±48 




HD 153948 


2453569.768 


1005 


-2 


-222 ± 82 


DD 




2453570.774 


942 


-2 


-185 ±89 


DD 




2453571.778 


901 


-3 


199 ± 90 




HD 317857 


2453571.810 


568 


-10 


-919 ± 20 


DD 


HD 318100 


2453569.804 


597 


-8 


-644 ± 59 


DD 


HD 162576 


2453568.870 


1861 


-14 


10 ±25 


DD 




2453569.828 


1568 


-13 


-19 ±28 






2453570.797 


2498 


-13 


13 ± 18 




HD 162588 


2453568.884 


1625 


-16 


93 ± 84 






2453569.846 


1039 


-15 


41 ± 141 






2453571.858 


2016 


-13 


36 ±70 




HD 162630 


2453570.819 


1229 


-15 


-61 ± 117 






2453571.872 


1407 


-13 


-157 ± 116 




HD 162656 


2453571.885 


1338 


+3 


24 ±78 




HD 162725 


2453566.876 


1654 


-14 


-101 ± 14 


DD 




2453569.865 


1608 


-13 


-21 ± 20 


DD 


HD 320764 


2453568.858 


1212 




-399 ± 408 






2453571.832 


1131 




137 ± 448 




HD 169842 


2453566.902 


1395 


-30 


180 ± 20 


DD 




2453569.890 


930 


-30 


230 ± 30 


DD 




2453570.858 


1441 


-30 


-153 ±21 


DD 


HD 169959 A 


2453570.885 


1378 


-53 


-691 ± 113 


DD 


HD 170054 


2453569.925 


781 


-28 


12 ± 108 






2453570.888 


1556 


-26 


-87 ± 54 




HD 170860 


2453571.901 


1005 


-13 


462 ± 236 


DD 


HD 172271 


2453567.897 


1441 


-21 


-265 ± 90 




HD 205073 


2453567.914 


1395 


-7 


4 ± 15 






2453569.972 


1102 


+ 18 


-8 ±21 




HD 205331 


2453567.922 


1614 


-5 


-49 ± 41 






2453570.003 


2022 


-4 


-9 ±33 






2453570.967 


1688 


-3 


143 ± 41 




BD +49-3789 


2453570.941 


528 


-6 


561 ± 400 




HIP 109911 


2453570.033 


545 


-10 


494 ± 230 
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Table 4. Comparison of ESPaDOnS and FORS 1 measurements of the same stars 



Star 


m v 


Teff 


v e sin i 


Instrument 


HJD 


At 


SNR 


(B z ) ± o-b 


fid 






(K) 


(kms" 1 ) 






(min) 




(G) 




HD 108945 


5.46 


8800 


65 


FORS1 




15 


1250 


-44+ 75 


nnn 










ESPaDOnS 


2453569.732 


12 


1880 


-111+ 28 


DD 










ESPaDOnS 


2453569.740 


12 


1870 


-75+ 29 


DD 










ESPaDOnS 


2453570.735 


15 


2490 


109 ± 20 


DD 


HD 153948 


9.35 


10600 


80 


FORS1 




70 


3100 


195 ± 53 


DnD 










ESPaDOnS 


2453569.768 


40 


1000 


-222+ 82 


DD 










ESPaDOnS 


2453570.774 


38 


940 


-185+ 89 


DD 










ESPaDOnS 


2453571.778 


37 


900 


199 ± 90 


ND 


HD 317857 


10.30 


10000: 


26 


FORS1 




75 


2150 


-1558± 54 


DDD 










ESPaDOnS 


2453571.810 


46 


570 


-919+ 20 


DD 


HD 318100 


9.84 


10600 


42 


FORS1 




47 


2100 


345 ± 64 


DnD 










ESPaDOnS 


2453569.804 


50 


600 


-644± 59 


DD 


HD 162725 


6.42 


9600 


32 


FORS1 




46 


4450 


-60± 33 


nnn 










ESPaDOnS 


2453566.876 


20 


1650 


-101+ 14 


DD 










ESPaDOnS 


2453569.865 


26 


1610 


-21+20 


DD 


HD 320764 


8.90 


8500 


225 


FORS1 




36 


2400 


-69+ 66 


nnn 










ESPaDOnS 


2453568.858 


33 


1212 


-399± 408 


ND 










ESPaDOnS 


2453571.832 


38 


1131 


137 ±448 


ND 


HD 169959A 


7.57 


11000 


57 


FORS1 




30 


2550 


-486+ 52 


DnD 










ESPaDOnS 


2453570.885 


25 


1378 


-691± 113 


DD 


HD 170054 


8.18 


14500 


26 


FORS1 




42 


2300 


64 ±68 


nnn 










FORS1 




36 


2200 


84 ±77 


nnn 










ESPaDOnS 


2453569.925 


52 


780 


12 ± 108 


ND 










ESPaDOnS 


2453570.888 


27 


1560 


-87+ 54 


ND 


HD 170860 


9.39 


13700 


85 


FORS1 




45 


2050 


-41+ 61 


nnn 










ESPaDOnS 


2453571.901 


43 


1005 


462 ± 236 


DD 



Table 5. Physical properties of stars that are probable or certain open cluster members, probable or certain Ap stars, and have 
ESPaDOnS magnetic field measurements. 



Cluster logf true DM Star logr c log L/L Q M/M Q fractional age S rms n B 









(yr) 






(K) 














(G) 




NGC 1039 


8.33 


+ 


0.06 


8.49 


HD 16605 


4.025 


1.65 


2.55 


± 0. 


15 


0.40 


±0.09 


1894 


5 


Coma Ber 


8.70 


± 


0.10 


4.70 


HD 108945 


3.944 


1.60 


2.30 


±0. 


15 


0.71 


±0.21 


93 


4 


NGC 6281 


8.45 


± 


0.10 


8.86 


HD 153948 


4.025 


1.86 


2.70 


±0. 


15 


0.62 


±0.18 


201 


4 


NGC 6405 


7.80 


± 


0.15 


8.44 


HD 318100 


4.025 


1.64 


2.50 


±0. 


10 


0.11 


±0.04 


517 


2 


NGC 6475 


8.41 


± 


0.05 


7.24 


HD 162576 


4.004 


2.17 


3.10 


±0. 


15 


0.81 


±0.17 


15 


3 












HD 162588 


4.033 


2.10 


3.05 


±0. 


20 


0.77 


±0.17 


34 


3 












HD 162725 


3.982 


2.36 


3.30 


± 0. 


20 


0.94 


±0.20 


69 


3 


NGC 6633 


8.75 


± 


0.05 


7.93 


HD 169842 


3.924 


1.62 


2.35 


± 0. 


10 


0.85 


±0.18 


315 


8 


IC 4725 


8.02 


± 


0.05 


10.44 


HD 170860 


4.137 


2.63 


4.25 


± 0. 


20 


0.71 


±0.15 


40 


1 


IC 4756 


8.44 


± 


0.06 


7.59 


HD 172271 


4.013 


1.98 


2.85 


± 0. 


15 


0.69 


±0.16 


265 


1 


NGC 7243 


8.06 


± 


0.10 


9.54 


BD+49 3789 


4.111 


2.23 


3.55 


±0. 


15 


0.51 


±0.15 


561 


1 












HIP 109911 


4.114 


2.29 


3.65 


±0. 


15 


0.54 


±0.16 


494 


1 



